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ABSTRACT 


Leaf-area index (LAI), the total one-sided surface area of leaf per ground surface area, is a 
key component of land surface models. We investigate the influence of differing, plausible 
LAI prescriptions on heat, moisture, and carbon fluxes simulated by the Community Atmo- 
sphere Biosphere Land Exchange (CABLEvl.4b) model over the Australian continent. A 
15-member ensemble monthly LAI data-set is generated using the MODIS LAI product and 
gridded observations of temperature and precipitation. Offline simulations lasting 29 years 
(1980-2008) are carried out at 25 km resolution with the composite monthly means from 
the MODIS LAI product (control simulation) and compared with simulations using each of 
the 15-member ensemble monthly- varying LAI data-sets generated. The imposed changes in 
LAI did not strongly influence the sensible and latent fluxes but the carbon fluxes were more 
strongly affected. Croplands showed the largest sensitivity in gross primary production with 
differences ranging from -90 to 60 %. PFTs with high absolute LAI and low inter- annual 
variability, such as evergreen broadleaf trees, showed the least response to the different LAI 
prescriptions, whilst those with lower absolute LAI and higher inter-annual variability, such 
as croplands, were more sensitive. We show that reliance on a single LAI prescription may 
not accurately reflect the uncertainty in the simulation of the terrestrial carbon fluxes, es- 
pecially for PFTs with high inter-annual variability. Our study highlights that the accurate 
representation of LAI in land surface models is key to the simulation of the terrestrial carbon 
cycle. Hence this will become critical in quantifying the uncertainty in future changes in 
primary production. 
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1. Introduction 


Land Surface Models (LSMs) describe the exchange of heat, moisture, and carbon be- 
tween the land surface and the atmosphere. There are a wide variety of LSMs used in both 
regional and global climate models, and they can vary considerably in complexity (Pitman 
2003). One key aspect which differentiates LSMs is whether they include phenology, and if 
dynamic, whether it is prescribed or simulated. In most LSMs, phenology is represented by 
the leaf area index (LAI), the total one-sided surface area of leaf per ground surface area. 

LAI is critical in any LSM since it affects the albedo of the terrestrial surface, and hence, 
the amount of net radiation available to drive sensible and latent heat. LAI also affects 
the partitioning of net radiation between sensible and latent heat fluxes (Verstraete and 
Dickinson 1986) because it controls the surface area of vegetation in direct contact with the 
atmosphere and affects the efficiency by which water can be transferred from within the 
vegetation to the atmosphere. Similarly, LAI affects the terrestrial carbon balance since it 
affects the photosynthesis and net primary productivity of a canopy. Finally, LAI influences 
rainfall interception and thereby affects the partitioning of rainfall between evaporation, 
throughfall, and runoff. 

The implementation of LAI in LSMs within regional and global climate models varies 
widely. At one end of the spectrum, some LSMs are coupled to dynamic vegetation models 
(e.g., Bonan et al. 2003), whereby LAI is a prognostic variable and responds to surface climate 
variations. However, climate biases from the regional and global atmospheric models make 
the realistic simulation of LAI difficult (Liu et al. 2008). As a consequence, most LSMs do 
not include dynamic vegetation and instead prescribe LAI. 

LAI can be prescribed according to plant functional types (PFTs) from look-up tables. 
These values are usually based on field observations and either held constant in time or 
allowed to vary seasonally. This method does not allow for inter-annual variability or vari- 
ations within PFTs; the same PFTs at different latitudes use the same LAI. Since this is 
not realistic, several studies have investigated the use of satellite derived LAI and shown 
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improvements in the simulation of surface climatology (e.g., Pielke et al. 1997; Buermann 
et al. 2001). The main impediment to the use of satellite derived LAI is the limited tempo- 
ral availability of these data. There is also an inherent assumption of stationarity for future 
climate simulations; the assumption that the present spatial and seasonal variations in LAI 
are representative of the future, even though they are clearly climate-dependant. 

Since LAI interacts with radiation, water balance and carbon balance it is a key parameter 
connecting the core components of climate and ecological modeling (Parton et al. 1996). One 
of the key characteristics of LAI is how it varies spatially (Bonan et al. 1993) and temporally. 
While LAI affects the interactions between the atmosphere at a point or grid scale (Bonan 
et al. 1993) this scales up to continental scales (Pitman et al. 1999) in uncoupled simulations. 
There is additional evidence that LAI affects the atmosphere at larger scales (Chase et al. 
1996). Most recently, van den Hurk et al. (2003) demonstrated that using remotely sensed 
LAI in a weather forecasting system affected the surface evaporation when evaporation 
formed a large term in the surface energy balance. They concluded that improved estimates 
of LAI could be an important method for improving model estimates of evaporation. 

The relationship between LAI and the terrestrial carbon balance is well documented from 
observational studies. Barr et al. (2004) investigated the influence of LAI on net ecosystem 
production in a deciduous forest in Canada and found a tight coupling between the annual 
maximum LAI and production. Saigusa et al. (2008) used data from flux-towers and found 
that temperate deciduous forests showed the greatest positive net ecosystem production after 
leaf expansion (higher LAI) in early summer. Duursma et al. (2009) used measurements from 
coniferous stands in Europe and found that LAI was a significant influence on gross primary 
production (GPP). Finally, Keith et al. (2012) used measurements at a single flux-tower 
site in Australia and focused on the carbon budget during drought years. They found that 
reductions in LAI due to insect attacks, in addition to drought stresses, contributed to a 
26% reduction in GPP and 9% reduction in ecosystem respiration as compared to years with 
drought stresses alone. 
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Some modelling studies have investigated the influences of vegetation parameters on the 
simulation of the terrestrial carbon fluxes and season length (e.g., White and Nemani 2003; 
Piao et al. 2007), but few explicitly focus on the influence of LAI versus meteorological 
forcing. This was recently investigated by Puma et al. (2013) in an offline LSM at four 
North American sites. They found that variations in LAI had a dominant control on GPP, a 
smaller but comparable effect on transpiration, a weak influence on total evapotranspiration, 
and a negligible impact on runoff. Additionally, they found that the effect of LAI on GPP 
is greater in energy-limited regimes as compared to moisture-limited regimes, except when 
vegetation exhibits little inter- annual variations in LAI. Hence, they conclude that an accu- 
rate representation of LAI inter-annual variability in LSMs is critical to accurately simulate 
GPP. 

Overall, it is clear that the way a LSM treats LAI is central to accurately simulating the 
heat, moisture, and carbon fluxes at the land surface. This paper focuses on the Community 
Atmosphere Biosphere Land Exchange Model (CABLE) (Wang et al. 2011). CABLE does 
not include a dynamic vegetation model by default, and hence the spatial and temporal 
variation of LAI are prescribed (prognostic LAI is implemented in later versions but not 
currently widely used). While several studies have used CABLE to answer wide-ranging 
research questions (e.g., Abramowitz and Gupta 2008; Cruz et al. 2010; Zhang et al. 2011b; 
Pitman et al. 2011; Wang et al. 2012; Exbrayat et al. 2012), only few studies have examined 
the influence of LAI on heat, moisture and carbon fluxes in CABLE. 

Zhang et al. (2013) ran global offline simulations with CABLE and conducted a sensitivity 
analysis by varying several vegetation and soil parameters, including LAI, by ± 50, 30, and 
20 % of the default values. Comparison of their simulations with other models (Rodell 
et al. 2004; Dirmeyer et al. 2006; Jung et al. 2009) showed that the influence of LAI was 
most noticeable in the middle and high latitudes of the northern hemisphere where broadleaf 
forests are the dominant PFT. However, Zhang et al. (2013) also point out that their imposed 
LAI perturbation do not necessarily reflect realistic uncertainties in estimates of LAI, and 
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additionally, only focussed on evapotranspiration and run-off. 

Lu et al. (2013) conducted an extensive parameter sensitivity analysis of CABLE over 
a single year at the global scale. They found that the at the global scale, the most impor- 
tant parameter affecting GPP is the maximum carboxylation rate, followed by LAI. When 
analysing each PFT separately, they also found LAI to be the second most important pa- 
rameter influencing GPP, except for evergreen broadleaf forests, whereby the initial slope of 
the response curve of potential electron was the second most important factor, followed by 
LAI. They carried out a similar analysis for latent heat, and found LAI to be the third most 
important factor globally, but results varied for each PFT. Namely, LAI was the most im- 
portant for deciduous needleleaf trees, second most important for evergreen needleleaf trees, 
third most important for evergreen broadleaf trees, deciduous broadleaf trees, and deciduous 
needleleaf trees, fourth most important for crops, and fifth most important for shrublands. 

Whilst the work of Zhang et al. (2013) and Lu et al. (2013) provide valuable insight into 
the sensitivity of CABLE to LAI, and it’s importance relative to other model parameters, 
the influence of realistic inter-annual variations in LAI on the surface energy and carbon 
balance remains un-known. This study provides a method of generating LAI ensembles, 
based on the MODIS LAI and the observed climatology, to address this knowledge gap. The 
next section describes the model set-up and the generation of the LAI ensemble. This is 
followed by an analysis of the influence of different monthly-varying LAI prescriptions on 
CABLE simulated surface energy and carbon fluxes. 


2. Methods 

a. Model Description 

CABLE is a LSM designed to simulate fluxes of energy, water and carbon at the land 
surface and can be run as an offline-model with prescribed meteorology (e.g., Wang et al. 
2011) or fully coupled to an atmospheric model within a global or regional climate model 
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(e.g., Mao et al. 2011). CABLE is a key part of the Australian Community Climate Earth 
System Simulator (ACCESS, see http://www.accessimulator.org.au), a fully coupled 
earth system science model, currently being used as part of the fifth assessment report of 
the International Panel on Climate Change. The version used in this study is CABLEvl.4b. 

In CABLEvl.4b, the one-layered, two-leaf canopy radiation module of Wang and Leuning 
(1998) is used for sunlit and shaded leaves and the canopy micro- meteorology module of 
Raupach (1994) is used for computing surface roughness length, zero-plane displacement 
height, and aerodynamic resistance. The model also consists of a surface flux module to 
compute the sensible and latent heat flux from the canopy and soil, the ground heat flux, as 
well as net photosynthesis. A soil module is used for the transfer of heat and water within 
the soil and snow, and an ecosystem carbon module based on Dickinson et al. (1998) is used 
for the terrestrial carbon cycle. A detailed description of each of the modules can be found 
in Kowalczyk et al. (2006) and Wang et al. (2011). 

LAI in CABLE is used to compute the roughness length of vegetation and the standard 
deviation of vertical velocities, which are used for the formulation of the aerodynamic resis- 
tances, and hence influence surface energy balance calculations. It is also used to compute 
the total flux density of radiation for sunlit and shaded leaves within the plant canopy radi- 
ation transfer model. This influences simulations of photosynthesis, stomatal conductance, 
leaf temperature, and energy and carbon fluxes as CABLE performs separate calculations 
for sunlit versus shaded leaves (Kowalczyk et al. 2006). Finally, LAI is used in the ecosys- 
tem carbon module where it directly influences GPP and autotrophic respiration (AR). 
Heterotrophic respiration (HR) is not directly driven by LAI, but by soil moisture and tem- 
perature. 

b. Model set-up 

CABLEvl.4b was used within the NASA Land Information System (LIS-6.1) (Kumar 
et al. 2006, 2008), a flexible software platform designed as a land surface modelling and 
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hydrological data assimilation system. A grid resolution of 0.25 x 0.25 degrees was utilised, 
covering continental Australia. The model was forced with the Modern Era Retrospective- 
analysis for Research and Applications (MERRA) reanalysis (Rienecker et al. 2011) at 3- 
liourly intervals and integrated from 1980-2008 and initialised from a previous 30-year spin- 
up. The forcing variables included incoming long-wave and short-wave radiation, air tem- 
perature, specific humidity, surface pressure, wind speed and precipitation. The MERRA 
reanalysis was bias-corrected for precipitation using the Australian Bureau of Meteorology 
Australian Water Availability gridded precipitation dataset (Jones et al. 2009), following 
Decker et al. (2012). Monthly ambient carbon-dioxide concentrations were prescribed using 
measurements from Baring Head, New Zealand (Keeling et al. 2005). 

In CABLEvl.4b, the background snow-free and vegetation- free soil albedo has to be 
prescribed. We used the MODIS derived snow-free background soil albedo data from Hould- 
croft et al. (2009). Bare soil regions, as defined by the IGBP land-use classification map 
(used in CABLE), are assigned the mean albedo over the data period (October 2002 to 
December 2006), whilst partially vegetated pixels are assigned a soil albedo derived from a 
linear relationship between albedo and the Normalised Difference Vegetation Index (NDVI). 
A linear regression model is then used to estimate the background soil albedo corresponding 
to zero green LAI (Houldcroft et al. 2009). The IGBP land-use classification was used, and 
radiative properties, including the leaf transmittance and reflectance values in the visible, 
near infra-red, and thermal regions were prescribed for each vegetation type following Avila 
et al. (2012). These values were obtained by adjusting estimates from Dorman and Sellers 
(1989) until the simulated albedo from CABLE closely approximated the MODIS observed 
broadband albedo. 

c. Simulations 

When running CABLE at a single site, LAI can be prescribed from observations at the 
site (e.g., Abramowitz and Gupta 2008; Wang et al. 2011; Li et al. 2012). When running 
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CABLE over a grid domain, LAI values are by default taken from a literature-based estimate 
for each PFT, and are fixed in time (e.g., Zhang et al. 2011a) or vary seasonally (Avila 
et al. 2012). For IPCC AR5 global climate simulations, the MODIS LAI product is used in 
CABLE within the ACCESS global circulation model. Since the aim of this paper is better 
inform the sensitivity of CABLE to LAI, we use the same MODIS LAI product (Yuan et al. 
2011) for our control simulation (1980-2008). This is carried out by prescribing monthly 
mean climatological LAI at each grid cell, based on monthly averages over the period of 
availability of the MODIS LAI data (2000-2008). 

To investigate the influence of LAI, a 15-member monthly-varying (1980-2008) LAI en- 
semble was generated using the MODIS LAI and gridded observations of maximum (Tmax) 
and minimum (Tmin) temperatures and precipitation from the Bureau of Meteorology Aus- 
tralian Water Availability Project (BAWAP) (Jones et al. 2009). The goal of reconstructing 
the LAI was to explore the model response to reasonable estimates of LAI variability and 
therefore, an ensemble approach based on simple linear regression between the MODIS LAI 
and the BAWAP data was used. 

The 8-day MODIS LAI was spatially aggregated from its original 0.05 by 0.05 degree grid 
to the BWAP 0.25 by 0.25 degree grid, by weighting each 0.05 cell by the area, summing 
the twenty-five 0.05 degree grid cells within each 0.25 cell, and finally normalizing by the 
total area within the course grid cell. This simple method avoids introducing unnecessary 
complexities that arise when the LAI is interpolated using subgrid scale plant functional type 
distributions. The 8-day, 0.25 degree fields where finally averaged to the monthly means by 
weighting each 8-day period according to the number of days from that time-span that fell 
within a given month. 

The 15-ensemble members were generated by linearly regressing the anomalous (found 
by removing the mean annual cycle) monthly MODIS LAI against Tmax, Tmin, and pre- 
cipitation from BAWAP at each 0.25° grid cell. The regressions were performed using data 
from the period 2000-2008, as this period is coincident with availability of the MODIS LAI. 
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The regressions were first performed separately for each variable and subsequently using all 
three variables to isolate the influence of each of Tmax, Tmin, and precipitation. Due to 
the lag between precipitation and vegetation greenness metrics in Southeastern Australia 
(Decker et ah 2012) we use a centered 5-point, linear regression, although similar results are 
obtained when only three points are included. The different sets of spatially distributed 
regression coefficients were calculated by randomly removing 25% of the data from each of 
the 15 regressions. 

Data was withheld as the data training period (2000-2008) occurs during a long-term, 
large scale drought in Australia. Limiting the temporal data in each of the regressions allows 
for uncertainty due to the training period selection and creates a larger spread among the 
final ensemble members. The 15 ensemble estimates of anomalous LAI were created by 
applying each of these 15 different, spatially explicit regression coefficients for the period 
1980-2008. A random Gaussian noise component with the mean and standard deviation 
given by the mean and standard deviation of the regression errors from each fitting was 
added during the construction of the LAI estimates. The added noise ensures that the errors 
associated with the fitting propagates to the final estimates, increases the spread between 
each of the ensemble members, and is consistent with the assumption that errors in LAI 
follow a Gaussian distribution (McColl et al. 2011). Finally these estimates of the LAI 
anomalies (constructed using all three data sources) were added to the mean annual cycle of 
the MODIS LAI to create the final LAI ensemble members. The spatially averaged ensemble 
spread of the anomalous LAI, relative to (i.e. divided by) the spatially averaged ensemble 
mean anomaly was 19.1% for the median, 22.9 % for the mean, 0.1 % for the minimum, 
and 133.6 % for the maximum. Whilst this range of LAI is smaller as compared to the 
range of LAI imposed by other studies, it suits the purpose of testing the influence of a 
climatologically driven LAI ensemble which is the aim of this study. 

Figure 1 shows the relationship between the MODIS LAI and the mean of the 15 member 
ensemble LAI reconstructions using only precipitation (Figure la), Tmax (Figure lb), Tmin 
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(Figure lc), and the combination of all three (Figure Id). The root mean square errors 
(RMSE) of the single variable regressions are 0.190, 0.194, and 0.200 respectively, while 
using all three variables results in a slightly better fitting (with an RMSE of 0.188). Figure 
1 demonstrates that while precipitation, Trnax, and Tmin, can be used to reconstruct the 
LAI, the slope of the fittings are less than one (0.982, 0.981, and 0.980, respectively). The 
combination of the three (Figure Id) yield a slope of 0.987, which is statistically larger 
than the slopes of the regressions using a single variable but still less than one. Due to 
the slightly better agreement with the MODIS observations for the period 2000-2008, the 
LAI reconstructed using all three variables was used for the model simulations. Overall the 
mean of the ensemble members reconstructs the LAI variability for the period 2000-2008 
with R2 values typically 0.3-0. 6, with some individual ensemble members better matching 
the observed LAI variability over this period. 

15 simulations were performed over this period using these monthly- varying LAI recon- 
structions. We note here that several studies on the influence of LAI on surface climatology 
use time- varying versus fixed LAI (e.g., van den Hurk et al. 2003) or apply a fixed fac- 
tor (e.g., double or half LAI, (Parton et al. 1996)). Since it is well established that the 
seasonal variation of LAI is not negligible (e.g., over croplands), and the use of remotely 
sensed LAI in LSMs generally improves surface climatology (Pielke et al. 1997; Buermann 
et al. 2001), we focus here on one of the most widely adopted remotely-sensed LAI prod- 
ucts, MODIS, and examine the sensitivity of CABLE to a MODIS-derived monthly varying 
ensemble LAI product, which is representative of the climatology. In summary, both the 
control and experiments are run over the same time-period, except that the control simula- 
tion has no inter-annual variation in LAI while the ensemble members are designed to reflect 
the climatology. 
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d. Data analysis 


The heat, moisture, and carbon fluxes were analysed separately for each dominant PFT, 
defined as PFTs with coverage greater than 1% of land points as shown in Figure 2. This was 
to avoid compensating effects between PFTs, as these have distinct seasonal signals as well 
as absolute magnitudes. For example, croplands, being a human-managed PFT, have higher 
seasonal variability than native vegetation. Additionally, the dense forested areas (evergreen 
broadleaf trees), have the highest absolute LAI, while most of inland Australia is sparsely 
vegetated with open shrublands with lower absolute LAI. Since the imposed changes in LAI 
are on the monthly time-scale, we compute monthly means and standard deviations of the 
fluxes and plot time series of the difference between the control and ensemble mean, with the 
standard deviation used to provide a measure of spread. Since the variations in the imposed 
LAI vary with time (monthly) and reflect the inter-annual variability in climatology inherent 
in the BAWAP gridded precipitation and temperature data-set, we perform a time-series 
rather than seasonal analysis (e.g., mean summer fluxes over the whole period). Additionally, 
we compute zero-lag cross-correlations between LAI and the fluxes to better quantify the 
response to changes in LAI. 


3. Results 

Figure 3 shows a monthly time-series of (a) the absolute (control-ensemble mean), and 
(b), percentage difference ((absolute difference/control) x 100) in LAI, heat, moisture, and 
carbon fluxes for open shrublands between 1980-2008. The zero-lag cross correlations with 
LAI are summarised in Table 1. The difference in LAI for open shrublands varies approxi- 
mately between -0.2 to 0.1, which represents a percentage change of -90 to 30 %. As expected, 
increases in LAI lead to a increase in vegetation transpiration (EV) and an decrease in soil 
evaporation (ES) as shown by the strong positive cross-correlation between LAI and EV and 
negative correlation with ES (Table 1). Although the absolute changes in EV are smaller 
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than ES, when expressed as a percentage change, they are larger by a factor of ~ 2-3. This 
is expected as the amount of leaf respiration is a direct function of LAI, whereas LAI only 
acts to partially inhibit soil evaporation. 

The effects of LAI on the absolute changes in mean monthly sensible (Qh) and latent 
(Qle) heat fluxes are small (< 1 W m -2 ), with percentage changes between -4 to 6 % only, 
and the correlations with LAI are lower as compared to EV and ES. These small changes in 
Qh and Qle corresponded with equally small changes in net radiation and surface albedo (not 
shown). Overall surface albedo in CABLE is a function of the vegetation albedo, background 
snow-free soil albedo, and snow albedo. The area covered by open shrublands is not densely 
vegetated, and hence it is the background soil albedo which largely determines the overall 
surface albedo. Thus, the relatively small perturbation in LAI imposed did not alter the 
overall surface albedo to a large extent and hence, the partitioning between Qh and Qle was 
not generally affected. 

The changes in the terrestrial carbon fluxes, on the other hand, showed a much stronger 
response to LAI. A decrease in LAI led to a decrease in autotrophic respiration (AR), and in- 
crease in heterotrophic respiration (HR), with strong positive cross-correlation between LAI 
and AR and weaker negative correlation with HR (Table l).When expressed as a percentage 
change, the differences in AR were up to 3-4 times larger than HR. This was expected, since 
HR is driven by below-canopy and soil processes, whilst AR is a direct function of LAI. 
Similarly, GPP was strongly positively correlated with LAI (we note that by convention in 
CABLE, downwards fluxes (i.e., GPP) are negative, but shown as positive here to remain 
consistent with the literature), as it is also a direct function of LAI, with percentage dif- 
ferences between -40 to 20 % (the same order of magnitude as the percentage change in 
LAI). 

For croplands (Figure 4), the absolute change in LAI varies between -0.6 and 0.6, cor- 
responding to a percentage change of approximately -160 to 40 %. This is larger when 
compared to open shrublands and all the other PFTs. Croplands, being a human-managed 
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PFT, have the highest seasonal and inter-annual variation in LAI (~ 0.3-1. 8) as compared 
to open shrublands (~ 0.3-0. 5) and the other PFTs, and hence the strongest response to 
monthly changes in precipitation, Tmax, and Tmin, which were used to generate the en- 
semble. The absolute changes in the heat and evaporative fluxes are an order of magnitude 
higher as compared to open shrublands (Figure 3), and the corresponding percentage changes 
are about double. Although the absolute changes in Qli and Qle are larger as compared to 
open shrublands, this change on a monthly time-scale is relatively small (the large percentage 
changes in Qh of up to 600 % still represent a small absolute change). The small absolute 
LAI of croplands is such that even large percentage changes did not change the surface 
albedo to a large enough extent to significantly alter net radiation. The absolute changes in 
AR, HR, and GPP are also an order of magnitude larger as compared to open shrublands, 
and the percentage changes are comparable to the imposed change in LAI. 

The changes for the other PFTs (woody savannas, savannas, and grasslands) showed 
similar trends (not shown), most noticeable in the carbon, rather than the turbulent heat 
fluxes. Evergreen broadleaf trees (Figure 5) had the smallest percentage change in LAI, 
since they have the largest absolute LAI values, and low inter-annual variability (~ 2. 8-3. 4). 
Hence, this PFT had the smallest response in the carbon fluxes (-4 to 6 %), with lower cross- 
corrrelations to LAI as compared to the other PFTs (Table 1). Evergreen broadleaf trees 
also showed a small positive correlation to HR of 0.46 (Table 1), whilst all other PFTs had 
a negative correlation, showing that a dense canopy can enhance HR. Another noticeable 
result for Evergreen broadleaf trees was that soil evaporation had a larger response to LAI 
as compared to vegetation transpiration in both absolute and percentage terms. This was 
a counter-intuitive result, as dense forested canopies would be expected to have a larger 
response of vegetation evaporation to LAI as compared to soil evaporation. To further 
investigate this, we conduced two extra simulations with large perturbations to the control 
LAI of ± 50 %. 

Figure 6 shows the seasonal difference in LAI imposed between the two experiments 
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(+50% minus -50%) and the subsequent changes to vegetation and soil evaporation (we 
show contours rather than time-series as the imposed LAI for these simulations has no inter- 
annual variability). As expected, a doubling of LAI results in an overall increase in vegetation 
transpiration and decrease in soil evaporation. However, the decrease in soil evaporation is 
almost twice as large as in the increase in vegetation transpiration, especially along the 
east coast where most Evergreen broadleaf trees are found. This is further demonstrated in 
Figure 7, showing the fraction of vegetation transpiration as a function of evapotranspiration 
(vegetation + soil) for both experiments. Over a semi- arid continent, changes in LAI result 
in a stronger response of soil evaporation as compared to vegetation transpiration. 

Whilst there are clear differences in the month-to-month variation of the heat, moisture, 
and carbon fluxes, increases in one period may be cancelled by a decrease later on. Addition- 
ally, we have not considered any spatial patterns in the changes in LAI and carbon fluxes. 
This is illustrated in Figure 8, showing the gridded cumulative monthly mean difference in 
LAI on carbon fluxes (cumulative changes in LAI < 5 have been masked out to highlight 
the largest changes). Clearly, the largest changes in LAI and carbon fluxes are restricted 
to the southeastern, rather than southwestern croplands (see Figure 2). This is due to the 
imposed change in LAI being almost twice as high for the southeastern as compared to the 
southwestern croplands, as illustrated in Figures 9a and 9b respectively. The larger response 
to LAI in southeast is due to the larger inter-annual variation in precipitation in this region, 
which was used to generate the LAI ensemble. 


4. Discussion 

The literature clearly suggests that the prescription of LAI in LSMs has a strong influence 
on the surface heat, moisture, and carbon fluxes. Hence we conducted a series of experiments 
to examine the influence of LAI variability in CABLE, as it is a widely used LSM in the 
Australian climate community and this sensitivity has not been previously tested. 
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Our results show relatively small impacts on the partitioning of available energy into the 
sensible and latent heat fluxes. Other studies have found much larger impacts, however, these 
were confined to regions of much larger changes in LAI compared to the changes imposed in 
this study. For example, Pitman et al. (1999) found large changes in total evaporative fluxes, 
but these were confined to regions where the absolute change in LAI was up to 3. Similarly, 
Bonan et al. (1993) found that LAI had a strong influence on the surface energy balance, but 
focussed on western US Conifer forests, the LAI of which varies from approximately 5 to 13. 
The imposed changes in LAI were much smaller in magnitude, but realistic and plausible 
, i.e., related to the climatology. Even when the LAI was doubled, the magnitude of the 
change was less than 1 for most of the continent (Fig. 6 (a)). Hence, the relatively small 
response of the evaporative fluxes is due to a small (but realistic) perturbation in LAI. 

The experiments with ± 50 % of the control LAI showed that doubling LAI resulted in a 
decrease in soil evaporation, which is twice as large as the increase in vegetation transpiration. 
This result is consistent with other studies which have shown that over half of the water lost 
through evapotranspiration over the Australian continent is through soil evaporation and by- 
passes plants almost entirely (Haverd et al. 2013). Similar results have been found elsewhere. 
Namely, van den Hurk et al. (2003) showed that in relatively dry (moisture limited) areas, 
where LAI values are relatively low, changes in LAI cannot result in large changes in surface 
heat and moisture fluxes as the land surface is already constrained by available soil water. 
In other words, variations in LAI cause the stronger response where surface evaporation uses 
a large proportion of the available energy. 

van den Hurk et al. (2003) did not allow for changes in LAI to alter the surface albedo, 
and hence, omitted a feedback important to our results. In our simulations, the variations 
in LAI imposed resulted in small changes in surface albedo, and subsequently small changes 
in net radiation. The small change in albedo is due to the relatively small perturbation in 
LAI imposed and because Australia is sparsely vegetated over large regions. It is therefore 
the background soil albedo, rather than the vegetation albedo, which has a large influence 
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on overall surface albedo in these regions. 

We found larger impacts on the terrestrial carbon balance, with LAI strongly positively 
correlated to GPP and AR, and negatively correlated with HR, consistent with both obser- 
vational (Barr et ah 2004; Saigusa et al. 2008; Duursma et ah 2009; Keith et ah 2012) and 
modelling (Puma et ah 2013) studies which report a tight coupling between LAI and primary 
production. This tight coupling is not unexpected as LAI is a key variable in the parameteri- 
sation of the carbon cycle. It determines not only the area of leaf that is potentially available 
to absorb light (and fix carbon via primary production, i.e., GPP), but also the amount of 
light attenuated and precipitation intercepted by the canopy. This in turn influences soil 
temperature, moisture, and evaporation, which drive heterotrophic respiration. However, of 
greater interest is the net ecosystem exchange (NEE) of carbon, i.e., the difference between 
GPP and the sum of HR and AR. If NEE in negative, then the land surface is a net source 
of carbon and a sink when positive. In all our simulations, NEE was always positive for both 
the control and the ensemble mean, and hence, the changes in LAI did not change the land 
surface to a source of carbon. 

The largest impacts were found for croplands, which have the highest inter-annual vari- 
ability in LAI. The changes were mostly restricted to the southeast, rather than southwest 
croplands, as the imposed changed in LAI was almost double in the former compared to the 
latter region. The southeast of Australia experiences higher inter-annual rainfall variability 
as compared to the southwest due to large-scale teleconnections (Risbey et al. 2009), and 
this signal was reflected in the LAI ensemble produced, as it is derived using gridded, sta- 
tion based precipitation and temperature data. The least impact was found for evergreen 
broadleaf trees, which had highest absolute LAI and lowest inter-annual variability. These 
results are consistent with Guillevic et al. (2002) and Puma et al. (2013), namely, that the 
impact of LAI variability is less for denser vegetation and moisture limited regions (low 
evaporative fraction). 

Whilst our results are broadly consistent with existing literature, they are constrained 
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by several caveats inherent of the study design. The model grid domain was restricted to 
Australia, due to the spatial extent of the BAWAP precipitation and temperature data used 
for generating the LAI ensemble, as well as bias correcting the forcing data. Hence our results 
are largely applicable to arid and/or semi-arid regions. Nonetheless, the results presented 
here should help inform the design of a broad range of future climate simulations whereby LAI 
is prescribed, especially when the focus is on the terrestrial carbon cycle. Our results are also 
limited to one particular LSM driven offline with a particular atmospheric forcing. Thus, our 
results would results would be worth extending via a multi-model evaluation of the sensitivity 
of LAI in LSMs that simulate the terrestrial carbon cycle. Despite inevitable caveats, our 
results highlight that the sensitivity testing of LSMs to LAI should be extended to include 
the terrestrial carbon cycle (rather than just heat and moisture fluxes). Additionally, the 
sensitivity of crop biomes to LAI highlights a need for the better representation of crop 
phenology in LSMs. This however remains a difficult challenge as crops, in contrast to other 
PFTs, are strongly and directly influenced by human intervention. 


5. Conclusions 

LAI is a critical component of any LSM. In this study, we performed a sensitivity anal- 
ysis of heat and carbon fluxes to perturbations in LAI using the CABLE LSM over the 
Australian continent on a monthly time-scale. We showed that whilst the influences of LAI 
perturbations on the heat and moisture fluxes were low, the impact on the terrestrial carbon 
balance was large, especially for croplands. Our results are consistent with earlier studies 
which have shown that PFTs with high inter-annual variability are the most sensitive to 
LAI perturbations, whilst dense vegetation is less sensitive, especially in moisture limited 
regimes. A key conclusion is therefore that care should be taken in accurately prescribing 
LAI, particularly when simulating the carbon cycle. Clearly, assigning fixed LAI to PFTs 
and/or using climatological means from remote sensing products, will not accurately reflect 
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the interannual variability of LAI which can have a large impact on the cumulative carbon 
fluxes. 

While our results focus on Australia, they provide several useful conclusions to the 
broader LSM community. First, using an ensemble of LAI products in simulations can 
be a very useful and straightforward method in establishing one element of uncertainty and 
the method used to generate the LAI ensemble here can be adapted to other regions and/or 
globally. Second, there is a clear need to assess the influence of LAI on the terrestrial carbon 
cycle at the global scale. To our knowledge, no studies have systematically addressed this 
issue, and this would provide a means to better quantify the uncertainty in future changes 
in the global terrestrial carbon cycle. Third, the sensitivities we find to LAI, particularly 
in respect of terrestrial carbon, points to the urgent need to resolve the parameterization 
of LAI more systematically in LSMs. Ideally, this is not through better prescriptions of 
LAI, rather it is via the addition of leaf phenology modules to LSMs. This highlights an 
important area of development in CABLE, as well as other LSMs which have no explicit 
dynamical representation of LAI. Finally, we also note that for a more complete assessment 
of the influence of LAI in LSMs, both the representation of vegetation through PFT maps 
and LAI variability should be analysed parallel to each other. 
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Table 1: Zero-lag cross-correlations between differences in leaf area index (LAI) and differ- 
ences in: vegetation transpiration (EV), soil evaporation (ES), sensible heat (Qh), latent 
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Figure 1: Scatter plot of the ensemble mean of the constructed LAI (m 2 m -2 ) versus the 
MODIS LAI (m 2 m -2 ) for each grid cell for the period 2000-2008 obtained using (a) precip- 
itation, (b) minimum temperature, (c) maximum temperature, and (d) precipitation, and 
minimum and maximum temperature. 
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2-Evergreen broadleaf (4.9 %) 
8-Woody savannas (5.2 %) 
10-Grasslands (3.4 %) 


7-Open shrublands (72.6 %) 
9-Savannas (4.8 %) 
12-Croplands (7.5%) 


Figure 2: Dominant plant functional types (PFTs), defined as greater than 1% of land points 
(masked inland regions in white are PFTs less than 1% of land points). 
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(a) Open shrublands (Absolute differences) 
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Figure 3: Time series of (a) monthly mean absolute differences, and (b) percentage differences 
(next page), in LAI, vegetation transpiration (EV), soil evaporation (ES), sensible heat (Qh), 
latent heat (Qle), autotrophic respiration (AR), heterotrophic respiration (HR), and gross 
primary production (GPP) between the control simulation and the ensemble mean for open 
shrublands (72.6 % of land points). The shaded region represents one standard deviation. 
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(b) Open shrublands (Percentage differences) 
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Figure 3: Continued 
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(a) Croplands (Absolute differences) 
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Figure 4: Same as in Figure 3 but for for croplands (7.5 % of land points). 
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(b) Croplands (Percentage differences) 



Figure 4: Continued 
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(a) Evergreen broadleaf trees (Absolute differences) 
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Figure 5: Same as in Figure 3 but for for evergreen broadleaf trees (4.9 % of land points). 
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(b) Evergreen broadleaf trees (Percentage differences) 
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Figure 5: Continued 
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Figure 6: Differences in (a) LAI, (b) vegetation evaporation (EV) (mm day" 1 ), and (c) 
soil evaporation (ES) (mm day" 1 ) between the experiment with +50% and -50% of the 
control LAI (the masked inland areas are regions where the gridded precipitation data used 
to generate the LAI ensemble was missing, and hence these points were excluded from all 
analysis for consistency). 
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Figure 7: Ratio of vegetation evaporation to total evapotranspiration (i.e., EV/(ES+EV)) 
for the experiments with (a) +50% of the control LAI, and (b) -50 % of the control LAI. 


40 


LAI 


AR 






-150 -100 -50 0 50 100 150 200 -1500 -1000 -500 0 500 1000 1500 


Figure 8: Gridded cumulative difference in monthly mean LAI and carbon fluxes 

(Gg month -1 ) between the control simulation and the ensemble mean (cumulative changes 
in LAI < 5 have been masked out to highlight the largest changes). 
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(a) Absolute differences (Croplands-SW) 
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Figure 9: Time series of monthly mean absolute differences in LAI, autotrophic respiration 
(AR), heterotrophic respiration (HR), and gross primary production (GPP) between the 
control simulation and the ensemble mean for (a) southwestern (SW) croplands, and (b) 
southeastern (SE) croplands (next page). 
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(b) Absolute differences (Croplands-SE) 
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Figure 9: Continued 
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